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In order to shed some light on the current discussion about /(i?)-gravity theories we derive and 
discuss the bounds imposed by the energy conditions on a general f(R) functional form. The null 
and strong energy conditions in this framework are derived from the Raychaudhuri's equation along 
with the requirement that gravity is attractive, whereas the weak and dominant energy conditions 
are stated from a comparison with the energy conditions that can be obtained in a direct approach 
via an effective energy-momentum tensor for /(i?)-gravity. As a concrete application of the energy 
conditions to locally homogeneous and isotropic /(i?)-cosmology, the recent estimated values of the 
deceleration and jerk parameters are used to examine the bounds from the weak energy condition 
on the parameters of two families of /(i?)-gravity theories. 

PACS numbers: 98.80.-k, 98.80.Jk, 04.20.-q 



I. INTRODUCTION 

The observed late-time acceleration of the Universe 
poses one of the greatest challenges theoretical physics 
has ever faced. In principle, this phenomenon may be 
the result of unknown physical processes involving cither 
modifications of gravitation theory or the existence of 
new fields in high energy physics. Although the latter 
route is most commonly used, following the former, an 
attractive and complementary approach to this problem, 
known as /(i?)-gravity [l[, examines the possibility of 
modifying Einstein's general relativity (GR) by adding 
terms proportional to powers of the Ricci scalar R to the 
Einstein-Hilbert Lagrangian 1 . 

Recently, several different functional forms for f(R) 
have been suggested in the literature (see, e.g., Ref. [3(]). 
These different /(i?)-gravity theories have also been dis- 
cussed in different contexts as, e.g., in the issues related 
to the stability conditions Q , inflationary epoch [B| , com- 
patibility with solar-system tests and galactic data Q, 
the late-time cosmological evolution [7(, among others. 
In the cosmological context, although these theories pro- 
vide an alternative way to explain the cosmic speed up 
without dark energy, the freedom in building different 
functional forms of f(R) gives rise to the problem of how 
to constrain from theoretical and/or observational as- 
pects these many possible /(i?)-gravities. Recently, this 
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1 Another interesting approach is related to the possible existence 
of extra dimensions, an idea that links cosmic acceleration with 
the hierarchy problem in high energyphysics, and gives rise to 
the so-called brane- world cosmology [2(. 



possibility has been explored by testing the cosmological 
viability of some specific forms of f(R) (see, e.g., Q). 

Additional constraints to f(R) theories may also come 
by imposing the so-called energy conditions [§, Q (for 
a pedagogical review, see Ref. [IP]). As is well known, 
these conditions were used in different contexts to de- 
rive general results that hold for a variety of situations. 
Thus, for example, the Hawking-Penrose singularity the- 
orems invoke the weak (WEC) and strong (SEC) energy 
conditions, whereas the proof of the second law of black 
hole thermodynamics requires the null energy condition 
(NEC). More recently, several authors have used the clas- 
sical energy conditions of GR to investigate some cosmo- 
logical issues, such as the phantom fields potentials 11 1, 
expansion history of the universe (l2| . as well as evo- 
lution of the deceleration parameter and their confront 
with supernovae observations [l3j ]. 

An important aspect worth emphasizing is that the en- 
ergy conditions were initially formulated in the context 
of GR 14]. In other words, this amounts to saying that 



one has to be cautious when using them in the more gen- 
eral framework, such as the /(i?)-gravity. On the other 
hand, we note that /(i?)-gravity theories share an inter- 
esting property: starting from the Jordan frame, where 
gravity is described only by the metric tensor, and mak- 
ing a suitable conformal transformation on the metric, 
it can be shown that any f(R) theory is mathematically 
equivalent to Einstein's gravity with a minimally coupled 
scalar field [lj|. 

Thus, in principle, one could think of "translating" the 
energy conditions directly from GR and impose them 
on the new effective pressure and effective energy den- 
sity defined in the Jordan frame. To test the viability 
of such a procedure, in this paper we derive the energy 
conditions for general /(ii)— gravity by using the physi- 
cal ultimate origin of the NEC and SEC, which is the 
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Raychaudhuri equation along with the requirement that 
gravity is attractive. We find that the NEC and the 
SEC for f(R), although similar, are different from that 
of Einstein's gravity. The resulting inequalities are then 
compared with what would be obtained by translating 
these energy conditions directly from the effective energy- 
momentum tensor for /(i?)-gravity in the Jordan frame. 
There emerges from this comparison a natural statement 
of the weak and dominant energy conditions in the con- 
text of the /(i?)-gravity theories. As a concrete appli- 
cation of the energy conditions for locally homogeneous 
and isotropic /(i?)-cosmology, we use the recent esti- 
mated values of the deceleration and jerk parameters, to 
examine the bounds from the WEC on the parameters of 
two families of f(R) -gravity. 



II. ENERGY CONDITIONS 

A. General Relativity 

Much of the technique we will use to derive the energy 
conditions for f(R) theories can be borrowed from the 
approach to the NEC and SEC in the context of GR. 
Therefore, for the sake of completeness, we shall first 
briefly review the approach to these energy conditions in 
GR. The ultimate origin of these energy conditions is the 
Raychaudhuri equation along with the requirement that 
gravity is attractive. To make clear this point, let u 11 
be the tangent vector field to a congruence of timelike 
geodesies in a spacetime manifold endowed with a metric 
g^v . Raychaudhuri's equation (see Ref. [16] for recent 
reviews) then reads 

^ = - 3 2 - + - R^u v , (1) 

where R^ is the Ricci tensor, and 9 , a 111 ' and u)^ v are, 
respectively, the expansion, shear and rotation associ- 
ated to the congruence defined by the vector field u M . 
This equation is a purely geometric statement, and as 
such it makes no reference to any gravitational field equa- 
tions. However, since the GR field equations relate R^ v 
to the energy- momentum tensor T^ v , the combination of 
Einstein's and Raychaudhy's equations can be used to 
restrict energy-momentum tensors on physical grounds. 
Indeed, since the shear is a "spatial" tensor one has 
a 2 = a^a^ v > 0, thus from Raychaudhury's equation 
it is clear for any hypersurface orthogonal congruences 
= 0) the condition for attractive gravity (conver- 
gence of timelike geodesies or geodesic focusing) reduces 
to R^u^u" > 0, which by virtue of Einstein's equation 
implies 

R^ u»u v = ( T„„ - ^g^J u»u» > , (2) 



where is the energy-momentum tensor, T is its trace, 
and where we have used units such that 8ttG = c = 1, 
which we shall adopt hereafter. Equation ([2|) is noth- 
ing but the SEC stated in a coordinate-invariant way in 
terms of T^ v and vector fields of fixed (timelike) charac- 
ter . In this way, the SEC in the GR context encodes 
the fact that gravity is attractive. In particular, note 
that for a perfect fluid of density p and pressure p 

T^ v = {p + p)u^u v -pg^ v , (3) 

and the restriction given by Eq. © takes the familiar 
form for the SEC, i.e, p + 2>p > 0. 

The evolution equation for the expansion of a congru- 
ence of null geodesies defined by a vector field fc M has 
the same form of the Raychaudhury equation, but with 
a factor 1/2 rather than 1/3, and with —R fJiV k fl k 1 ' as the 
last term (see Ref. [1C§ for more details). In this case, 
the condition for the convergence (geodesic focusing) of 
hypersurface orthogonal (uj^Z = 0) congruences of null 
geodesies along with Einsteins's equation implies 

R^k v = T^k v >0, (4) 

which is the NEC written in a coordinate-invariant 
way Thus, in GR the NEC ultimately encodes the 
null geodesic focusing due to the gravitational attraction. 
We note that for the energy-momentum tensor of a per- 
fect fluid [Eq. Q] it reduces to well-known form of the 
NEC, i.e., p + p > 0. 

B. /(i?.)-Gravity 

Since the Raychaudhuri equation and its null-geodesics 
counterpart hold for any geometrical theory of gravita- 
tion, we will keep the above physical motivation for the 
SEC and NEC in the /(i?)-gravity context. To this end, 
we shall use the above Raychaudhuri's geometric rela- 
tions [see Eqs. ([2]) and ([4])] along with the attractive char- 
acter of gravitational interaction to derive these energy 
conditions in the context of /(i?)-gravity theories. 

We begin by recalling that the action that defines an 
/(i?)-gravity is given by 

S = J d A x^gf(R) + S m (5) 

where g is the trace of the metric g^ u , R is the Ricci 
scalar and S m the standard action for the matter fields. 
Varying this action with respect to the metric we obtain 
the field equations 

f'R»u - ~g» v - (v M v„ - 9flv □ ) /' = Tftv , (6) 

where a prime denotes differentiation with respect to R 
and □ = g a & Vq.V/3 . In order to use an approach to the 
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NEC and SEC similar to that in GR context, we note 
that Eq. © can be rewritten as 

T 

R^v = Tfj,iy 7^9liv 7 (7) 

where i 

T = jj(T + f -Rf- 3D f) , 

V = j t [T^ + (V M V„ - ^D)/'] , 

and T = g^T^ . 

Now, for the homogeneous and isotropic Friedmann- 
Lemaitre-Robertson- Walker (FLRW) metric with scale 
factor a(t), the relations R^k^k" > and R^ v u^u v > 
[see Eqs. j2]) and ([2])] along with Eq. ([7]) lead, for a perfect 
fluid T^ v , to the following inequalities: 

p+p>0, and (R- RH) f" + R 2 f" >0, (8) 



P + 3 P - f + Rf + 3[R + RH) f" + 3R 2 f" > , (9) 



where a dot denotes derivative with respect to time, 
H = a/a is the Hubble parameter, and we have made 
the assumption that /' > 0, which ensures the regularity 
of the conformal rescaling [l5(. Thus, the inequalities ([8]) 
and © are, respectively, the NEC and SEC requirements 
in the context of /(i?)-gravity for a perfect fluid FLRW 
model. We note that the well-known forms for the NEC 
(p + P > 0) and SEC (p + 3p > 0) in the context GR 
can be recovered as a particular case of these energy con- 
ditions in f(R) theories for / = R, as one would have 
expected. 

An important point to derive the week and domi- 
nant energy conditions (WEC and DEC, respectively) 
in /(i?)-gravity is the fact that the above NEC and SEC 
[Eqs. @ and ([9])] can also be recovered as an extension of 
these conditions in GR. In fact, in /(i?)-gravity theories 
one can define an effective energy-momentum tensor as 
(see, e.g., Ref. [l7|) 



T eff = — 



Tfjtv + -(/ - Rf)g^ + (V M V„ - g^n)f 



(10) 

from which one defines an effective energy density and 
pressure by 



1 

T 



P+\(f-Rf)-3RHf 



(11) 



and 

Pes 
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P-^(f- R f) + ( R + 2RH)f" + R 2 f" 



(12) 

which in turn make apparent that the NEC and WEC 
given by Eqs. ([8]) and ([9]) can be obtained in a similar way 
as that in GR context. Thus, following and extending the 



GR approach to WEC and DEC, we have that in f(R)- 
gravity theories, in addition to Eq. ©, the WEC requires 
that 

P+l(f-Rf')-3RHf">0, (13) 

whereas the DEC fulfillment, besides the inequalities ([8]) 
and (Tl3l). reads 



p — p + f — Rf — (R+ 5RH)f" - R 2 f" > . (14) 



As one may easily check, for / = R, Eqs. ([13]) and (fT4[) 
give p > and p—p > 0, whose combination with Eqs ([5]) 
give, respectively, the well-known forms of the WEC and 
DEC in GR (see,e.g., Ref. [I3|)- We also note that ac- 
cording to Ref. [18j | if the energy momentum tensor is 
trace- free, then any homogeneous and isotropic solution 
of GR is also a particular solution of a /(i?)-theory pro- 
vided /(0) = , and /'(0) ^ 0. Thus, it is clear from the 
above Eqs. ©, ©, (13]) and ([14]) that in these cases the 
GR energy conditions are recovered. 

For the sake of completeness, we note that rather 
than the geodesic focusing the so-called sudden future 
singularity (SFS) is another type of singularity that 
arises in the context of perfect fluid FLRW cosmol- 
ogy 0, 0, HH, H3, HH . It comes about in the absence of 
an precise equation of state and is a singularity of pres- 
sure p only, with finite energy density p. In the context of 
GR, the SFS can occur with no violation of the SEC and 
WEC but with violation of the DEC [H, H3, HH, E3, EH . 
The nature of SFS in terms of geodesic completeness has 
been discussed in Ref. 24]. A discussion of the SFS in 
the context of the /(i?)-gravity, however, is beyond the 
scope of the present paper. In this regards, we refer the 
readers to Ref. [2(| (see also Ref. [25j|). 



III. CONSTRAINING f(R) THEORIES 

The energy-condition inequalities (O, ©, (fT^j) 
and (TTJJ can also be used to place bounds on a given 
f(R) in the context of FLRW models. To investigate 
such bounds, we first note that the Ricci scalar and its 
derivatives for a spatially flat FLRW geometry can be 
expressed in terms of the deceleration (g), jerk (j) and 
snap (s) parameters, i.e, 

R = -6H 2 {1 - q) , (15a) 
R = -6H 3 {j -q-2) , (15b) 
R = -QH 4 (s + q 2 + 8q + 6) , (15c) 



where I2i 



1 d 



3 



1 

1p7 



and 



1 "a 
1F~ 



(16) 



In terms of the present-day values for the above param- 
eters, Eqs. ©, ©, H3]) and (O can be rewritten as 
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(NEC) 



Po + Po > and 



[•so - Jo + (go + l)(ffl> + 8)]fo + G[Ho(jo - (to - 2)] 2 /o" > , 



(17a) 



(WEC) 



2 Po + /o + 6H 2 (1 - gb)/5 + 36ff 4 (io - <to - 2)/£ > , 



(17b) 



(SEC) po + 3 Po + /o - 6# 2 (1 - go)f{, - +6H*(s + j + q 2 + 7q + 4)/ ' + 3[6tf 3 (j -q - 2)] 2 / " > , (17c) 



(DEC) 



Po 



Po + 6H 2 (l - - 6tf 4 [s + (9o - l)(go + 4) + 5j ]/ 7 - [6H$(J -q - 2)] 2 / " > . (17d) 



r 



A. f(R) = R + aR n 

To exemplify how the above conditions can be used to 
place bounds on f(R) theories, we first note that, apart 
from the WEC [Eq. (jT7bJ)] , all above inequalities depend 
on the current value of the snap parameter sq- Therefore, 
since no reliable measurement of this parameter has been 
reported hitherto, in what follows we shall focus on the 
WEC requirement in the confrontation of the energy con- 
dition bounds in f(R) — gravity with observational data. 

As a first concrete example, we shall consider the fam- 
ily of theories with f(R) of the form 



f(R) = R + aR n 



(18) 



where n is an integer and a is a constant that can assume 
positive or negative values. For this class of theories we 
note that the constraints from WEC depend only on ob- 
servational values qo and jo > besides n and a that specify 
a particular theory. The inequality for WEC fulfillment 
condition (|17b[) can be written in terms of these param- 
eters as 



a{-l) n (An 2 - (A + l)n + 1) > . 



(19) 



where A = (jo — qo — 2)/(l — qo) 2 - In what follows 



we consider qn 



-0.81± 0.14 and j = 2.16 



+0.81 
-0.75' 



as 

given in Ref. [27j. The roots of the quadratic function 
are (1, l/A), so it has positive values for 1 > n > 3.377 
and negative values for 1 < n < 3.377. 

Inequation (TT5]) makes apparent that the observance 
of the WEC clearly depend upon the sign of a and the 
values of n. In what follows we consider the two possible 
signs for a along with the requirement that f'(R) > 
for all R, which reduces to 



-1)™ n(3.3# 



,2n-2 



< 1 



(20) 



for f(R) given by |IS]). Thus, for the family of f(R) 
theories given by Eq. (TTB")) the WEC is obeyed in the 
following cases: 

(i) a > . Here the allowed values for n 
can be grouped in following sets: n = 
{3,1,-2,-4,-6,...} and n = {4,6,8,...}, with 
< a < [n(3.3ff ) 2n_2 ] _1 for this last set. 



(ii) a < . In this case the allowed sets are n = 
{2, -1, -3, -5, . . .} and n = {1, 5, 7, 9, . . .}, with a 
within the range -[n(3.3i?o) 2 " _2 ] _1 < a < for 
this last set. . 

As a concrete application, it is clear from the above 
analysis that the so-called quadratic gravity, i.e. R + aR 2 
gravity (see, e.g., Ref. [§]) obey the WEC only if a is 
negative. Another interesting example is the / = R — 
/i 4 /i? gravity theory (see Ref. (28|), which correspond to 
WEC fulfillment case a < 0, n = — 1. The closely related 
/ = R — /i 6 / R 2 theories, however, violates the WEC, 
since it corresponds to a < and n — — 2. We also note 
that the violation of WEC by a specific member of a class 
of /(-R)-gravity theories clearly lead to the breakdown of 
the DEC. 



B. 



f(R) = aR n 



This type of /(i?)-gravity theories has been proposed 
in the investigation of galactic environments |29l . [30j ] . 
gravitational waves and lensing effects [3l|. as well as 
FLRW dynamics via phase space analysis [32J, [33 1 . For 



this class of theories the WEC constraints are again given 
by Eq. ([19]). but the requirement that f'(R) > (for all 
R) now takes the form 



«(-!)" n(3.3iJ ) 2 



< 0. 



(21) 



As before we group this class of gravity theories in two 
cases depending on the sign of a. We find that, for in- 
tegers values of n, the WEC is obeyed in the following 
cases: 

(i) a > . The allowed values for n are n = 
{3,1,-2,-4,-6,...}. 

(ii) a < 0. In this case we have the set n — 
{2,-1,-3,-5,...}. 

Finally, we note that the WEC- fulfillment case (i) may 
accommodate the observational value n — 3.5 obtained 
by fitting the Supernovae la Hubble diagram with this 
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type of power-law gravity [30| . In Ref. [11] the au- 
thors found that this type of theories admit simple exact 
solutions for FLRW models and that a stable matter- 
dominated period of evolution requires n > 1 or n < 3/4. 
Clearly, the first constraint do not violate the WEC, 
while the second do. 

IV. FINAL REMARKS 

/(-R)-gravity provides an alternative way to explain 
the current cosmic acceleration with no need of invoking 
either the existence of an extra spatial dimension or an 
exotic component of dark energy. However, the arbitrari- 
ness in the choice of different functional forms of f(R) 
gives rise to the problem of how to constrain the many 
possible /(i?)— gravity theories on physical grounds. 

In this paper we have shed some light on this question 
by discussing some constraints on general f(R)— gravity 
from the so-called energy conditions. Starting from the 
Raychaudhuri's equation along with the requirement that 
gravity is attractive, we have derived the null and strong 
energy conditions in the framework of /(ii)— gravity and 
shown that, although similar, they differ from those de- 
rived in the context of general relativity. By comparing 
the NEC and SEC inequalities [Eqs. © and ©] with 
what would be obtained by translating these energy con- 
ditions directly from the effective energy-momentum ten- 
sor for /(i?)-gravity in the Jordan frame, we also have 
obtained the general expressions for the weak and dom- 
inant energy conditions [Eqs. (Ti~3")) and l|14p]. As a con- 



crete example of how these energy conditions require- 
ments may constrain /(i?)-gravity theories, we have dis- 
cussed the WEC bounds on two different f(R) classes, 
i.e., f(R) = R + aR n and f(R) = aR n . An interesting 
outcome of this analysis is that the so-called quadratic 
gravity (R + aR 2 ), originally discussed in the context of 
primordial inflation, violates the WEC requirement for 
positive values of a . 

Finally, we emphasized that although the energy con- 
ditions in /(i?)-gravity discussed in this article has 
well-founded physical motivation (Raychadhuri's equa- 
tion along with the attractive character of gravity) the 
question as to whether they should be applied to any 
solution of /(i?)-gravity theories is an open question, 
which is ultimately related to the confrontation between 
theory and observations. We recall, however, that this 
confrontation in GR context indicate that all energy con- 
ditions seem to have been violated at some point of the 
recent past of cosmic evolution [l2| ■ 
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